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I.  INTRODUCTION 


.  , 

Structural  analysis  of  projectiles  and  other  ballistic 
components  often  require  the  modelling  of  screw  threads. 
Typically,  many  finite  elements  are  required  in  the  immediate 
vicinity  of  the  joint  to  accurately  define  and  model  a 
thread.  In  addition,  reconciling  the  effect  of  dimensional 
tolerances  of  threads  is  not  a  simple  effort.  For  example,  if 
an  analyst  were  required  to  perform  a  complete  analysis  of  a 
threaded  region,  including  tolerance  considerations,  at  least 
four  iterations  would  be  required.  These  would  include  the 
mating  of  large  and  small  pitch  diameter  male  threads  (within 
the  tolerance  range)  with  both  large  and  small  female  threads 
(also  within  tolerances).  .^2T7 

A  new  technique  for  modelling  threads  was  developed 
during  a  study  of  ogive  failures  of  the  M483  155-mm  artillery 
projectile  to  simulate  the  coupling  of  the  aluminum  ogive  to 
the  steel  body.  The  technique  employs  a  simple  matrix  of  gap 
elements  to  model  thread  behavior  as  opposed  to  the  labor- 
intensive  method  of  representing  threads  with  a  finite 
element  mesh.  Substantial  time  and  computer  memory  savings 
can  be  realized  by  employing  this  technique  with  no  penalty 
in  accuracy. 


These  new  "thread  elements"  greatly  simplified  the 
effort.  Otherwise,  it  would  have  been  necessary  to  use  very 
small  elements  to  accurately  grid  each  of  the  threads  in  the 
joint  region.  For  a  dynamic  analysis  in  which  stress  waves 
were  a  consideration,  this  would  have  required  a  very  small 
time  step  which  is  proportional  to  element  size  in  this  case. 
The  combination  of  small  time  steps  and  numerous  elements 
would  generate  output  data  that  may  exceed  the  limits  of  the 
computer  storage  device.  Also,  the  run  times  would  be 
prohibitively  long.  For  example,  it  requires  one  thread 
element  (set  of  three  gap  elements)  to  model  one  thread  pair 
and  at  least  eight  material  elements  to  realistically  mesh  a 
single  male  or  female  thread.  Also,  slide  lines  would  have  to 
be  included  at  the  closing  flanks  which  is  built  into  the 
thread  elements.  The  new  technique  is  also  useful  in  the 
study  of  thread  fit  and  tolerances.  The  thread  elements  can 
be  used  to  model  these  effects  on  the  behavior  of  the  overall 
structure  by  adjusting  the  "modulus"  in  the  axial  direction. 
This  modulus  might  be  determined  from  an  experimental  pull 
test  on  a  thread  specimen  where  a  stress-strain  curve  can  be 
developed  for  the  thread  "material".  The  modulus  determined 
from  such  a  test  can  then  be  input  directly  as  a  material 
property  for  the  gap  elements  used  to  model  threads.  A  set  of 
three  such  gap  elements  are  required  to  construct  a  standard 
thread  element.  The  approach  here  has  many  benefits  in 
analyzing  coupled  bodies.  However,  if  the  details  of  the 
stress  and  strain  state  in  the  vicinity  of  the  threads  must 
be  studied,  then  there  is  no  recourse  other  than  the 
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construction  of  a  detailed  local  model. 


II.  METHODOLOGY 

The  ANSYS^  finite  element  program  is  used  for  all 
analyses  included  in  this  report.  However,  any  general  finite 
element  code  with  gap  element  capability  could  be  used.  The 
user  must  know  a  priori  that  small  deformations  are  expected 
as  there  is  no  provision  for  releasing  axial  constraints 
should  the  threads  completely  decouple  radially.  The  threaded 
joint  used  to  develop  this  technique  is  found  on  the  M483 
155-mm  artillery  projectile  between  the  steel  body  and 
aluminum  ogive.  An  ongoing  study  of  ogive  failures  provided 
the  opportunity  for  the  development  of  this  technique.  The 
original  study  of  these  failures  was  performed  by  Southwest 
Research  Institute^  whose  finite  element  analysis  includes 
gridding  of  the  threaded  region.  Their  analysis  provides  the 
"geometry”  technique  as  a  basis  of  comparison  to  the  present 
"material  model"  technique  for  representing  thread  behavior. 

A  threaded  joint  has  basically  four  possible  responses 
to  any  given  load.  It  can  be  sheared  in  tension  or 
compression  in  response  to  a  tensile  or  compressive  axial 
load  respectively.  Also,  the  threads  can  either  compress 
together  or  separate  in  response  to  a  compressive  or  tensile 
radial  load.  The  joint  studied  in  this  case  has  a  buttress 
thread  since  the  expected  axial  load  can  only  act  in  one 
direction.  A  section  of  rotation  of  the  threaded  region  is 
illustrated  in  Figure  1.  There  is  virtually  no  chance  that 
non-load  bearing  flanks  will  meet  in  this  type  of  thread 
series. 

A  schematic  representation  of  a  buttress  thread  is  shown 
in  Figure  2  along  with  its  finite  element  model.  The  flank 
gap  element  angle,  o,  relative  to  the  axis  is  defined  to  be 
the  same  as  the  actual  flank  angle  and  there  is  one  radial 
gap  element  for  each  thread. 

A  gap  element  has  only  two  properties:  a  stiffness 
modulus  and  coefficient  of  friction.  The  modulus  of  the 
closed  gap  elements,  is  calculated  by  adding  reciprocal¬ 
ly  the  moduli  of  the  two  mating  materials  (the  same  as  adding 
two  spring  constants  in  series,  E^-  »  +  E2)).  Since 
each  "spring"  is  only  half  the  length  typfcally  used  in  this 
well-known  equation,  the  modulus  must  be  multiplied  by  2. 
Thus,  a  gap  element  between  aluminum  and  steel  would  have  a 


[2]  Cox,  P.A.,  et  al,  "A  Root  Cause  Analysis  of  the  M483A1 
Ogive,"  Southwest  Research  Institute,  April  1987. 
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Figure  2 .  Schematic  Representation  of  a  Buttress  Thread  and 
its  Finite  Element  Model. 
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modulus  of  15  Mpsi.  Sliding  is  permitted  after  closure,  and  a 
coefficient  of  friction  can  be  selected  from  a  handbook  or 
determined  from  other  means. 

Two  finite  element  models  of  the  threaded  region  are 
illustrated  for  comparison  in  Figure  3.  The  actual  threads 
are  modeled  in  the  grid  on  the  left  and  thread  elements  are 
depicted  in  the  grid  on  the  right.  The  dotted  line  on  the 
left  model  represents  the  surface  where  internal  pressure  of 
12,000  psi  is  applied. 


III.  RESULTS 

One  of  the  responses  of  the  joint  which  can  be  measured 
is  the  separation  of  the  closing  faces.  This  gives  an 
indication  of  the  overall  strength  of  the  threads.  The 
Southwest  Research  analysis  predicts  that  the  butt  joint 
between  the  body  and  ogive  will  open  .009  inch  at  peak  load 
while  the  thread  element  model  predicts  a  .010  inch 
separation  indicating  an  agreement  to  within  ten  percent. 

The  gap  separation  indicated  on  the  ANSYS  color  stress 
plot  in  Figure  4  compares  favorably  (within  ten  percent)  with 
the  plot  above  it  of  gap  separation  versus  time  from  SRI.  As 
far  as  stress  levels  in  thread  roots  and  crests  are  concerned 
the  thread  element  model  offers  no  information. 


IV.  DISCUSSION 

The  thread  element  model  represents  a  method  by  which 
the  analyst  can  more  easily  develop  the  finite  element  grid 
for  jointed  structures.  A  simple  model  is  more  easily 
understood  both  by  the  analyst  and  those  to  whom  the  results 
are  conveyed.  Also,  the  chances  of  modelling  errors  are 
reduced  with  less  complex  models. 

In  addition  to  model  simplicity,  less  time  is  required 
for  grid  development  and  computer  computations,  often  a 
concern  with  small  networks.  Also,  fewer  elements  mean  less 
computer  storage  space  is  needed.  In  some  cases,  the  time  and 
memory  reductions  can  be  significant.  This  is  best 
illustrated  when  analyzing  a  structure  for  stress  waves  or 
vibration^ . 


[3]  Bender,  J.M.  "Stress  Waves  Through  Joints  of  Artillery 
Projectiles,"  Ballistic  Research  Laboratory-TR-2980 ,  Jan  1989. 
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Figure  3.  Comparison  of  Finite  Element  Models. 


When  performing  a  finite  element  dynamic  analysis,  the 
time  step  is  a  function  of  the  element  size^.  The  element 
size  for  a  thread  element  is  the  distance  between  two 
threads.  In  contrast,  the  geometry  method  may  require  the  use 
of  many  elements  to  span  a  thread.  Consequently,  the  time 
step  must  be  much  smaller  to  accomodate  the  smaller  elements 
for  accurate  modelling  of  dynamic  effects  such  as  vibration 
and  stress  waves.  The  disadvantage  with  smaller  time  steps 
is  that  many  more  iterations  are  required  during  the  time 
span  of  the  analysis.  For  example,  if  element  size  is  reduced 
to  one  fourth  of  its  original  size,  the  time  step  is  likewise 
reduced.  This  means  it  will  take  four  times  the  number  of 
time  steps  to  span  the  duration  of  the  analysis  since  each 
time  step  represents  one  computation  of  strains.  Almost 
always,  this  means  a  similar  increase  in  computer  time  and 
memory  storage. 


[4]  Kohnke,  P.C.,  ANSYS  Theoretical  Manual ,  Swanson 
Analysis  Systems,  Inc.,  1983. 
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Washington,  DC  20360 

1  Commander 

Naval  Sea  Systems  Command 
Washington,  DC  20362 

1  Commander 

Naval  Ship  Systems  Command 
Washington,  DC  20360 

1  Commander 

Naval  Ship  Engr.  Center 
Hyattsville,  MD  20782 

1  Superintendent 

Naval  Postgraduate  School 
ATTN;  Dir  of  Lib 
Monterey,  CA  93940 

1  Commander 

Naval  Air  Development  Center 
Johnsville 

Warminster,  PA  18974 
1  Commander 

Pacific  Naval  Missile  Center 
Point  Mugu,  CA  95042 

1  Commander 

Naval  Research  Laboratory 
Washington,  DC  20375 


No.  Of 
Copies 


DISTRIBUTION  LIST 


Organization 

2  Commander 

Naval  Ship  Research  and 
Development  Center 
Bethesda,  MD  20084 

4  Commander 

Naval  Research  Laboratory 
ATTN :  W . J .  Ferguson 

C.  Sanday 
H.  Pusey 
F.  Rosenthal 
Washington,  DC  20375 

2  Commander 

Naval  Surface  Weapons  Center 
ATTN:  G-13,  W.D.  Ralph 

Dahlgren,  VA  22448 

1  Commander 

Naval  Surface  Weapons  Center 
ATTN:  Code  X211,  Lib 

Silver  Spring,  MD  20910 

1  Commander 

Naval  Surface  Weapons  Center 
ATTN:  Code  G-33,  T.N.  Tschirn 
Dahlgren,  VA  22448 

2  Commander 

Naval  Weapons  Center 
ATTN:  Code  3835,  R.  Sewell 
Code  3431,  Tech  Lib 
China  Lake,  CA  93555 

3^  Commander 

Naval  Ordnance  Station  Indian  Head 

ATTN:  Tech  Library 

Indian  Head,  MD  20640-5000 

2  Commandant 

USA  Infantry  School 
ATTN:  ATSH-CD-CSO-OR 

Fort  Banning,  GA  31905 

2  AFRPL 

ATTN:  W.  Andrepont 
T .  Park 

Edwards  AFB,  CA  93523 


No.  Of 

Copies  Organization 

1  AFOSR 

ATTN:  W.J.  Walker 

Bolling  AFB, 

Washington,  DC  20332 

3  AFATL  (DLA) 

Gun  and  Rocket  Division 
ATTN:  W.  Dittrich,  DLD 

D.  Davis,  DLDL 
Eglin  AFB,  FL  32542 

1  Air  Force  Armament  Lab 
ATTN:  AFTL  DLODL 
Eglin  AFB,  FL  32542 

1  Gun  Test  Branch 

AD  3246  Test  W/TETFG 
Eglin  AFB,  FL  32542 

2  AFWL/SUL 

Kirtland  AFB,  NM  87117 

2  ASD 

Wright-Patterson  AFB, 

OH  45433 

1  ASD  (ASAMCG)  Tech  Lib 
W-P  AFB,  OH  45433 

1  AFML  (T.  Nicholas) 

W-P  AFB,  OH  45433 

1  Director 

U.S.  Bureau  of  Mines 
Twin  Cities  Research  Center 
P.O.  Box  1660 
Minneapolis,  MN  55111 

2  Battel le  Columbus  Labs 
ATTN :  L . E .  Hulber t 

J.E.  Backofen,  Jr. 
505  King  Avenue 
Columbus,  OH  43201 

2  Battelle  Pacific  NW  Lab 
ATTN:  M.  Garnich 
M.  Smith 
P.O.  Box  999 
Richland,  WA  99352 
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DISTRIBUTION  LIST 
No.  of 

Copies  Organization 


No.  of 

Copies  Organization 

1  HQDA  (SARD-TR) 

Washington,  DC  20310-0001 

1  Commander 

US  Army  Laboratory  Command 
ATTN:  AMSLC-DL 
Adelphi,  MD  20783-1145 

1  Commander 

Armament  RD&E  Center 
US  Army  AMCCOM 
ATTN:  SMCAR-MSI 

Picatinny  Arsenal, 

NJ  07806-5000 

2  Commander 

Armament  RD&E  Center 
ATTN:  SMCAR-TDC 
Picatinny  Arsenal, 

NJ  07806-5000 

1  Commander 

USA  Aviation  Systems  Cmd 
ATTN:  T^SAV-DACL 
4300  Goodfellow  Blvd. 

St.  Louis,  MO  63120-1798 

1  Director 

US  Army  Aviation  Research 
and  Technology  Activity 
Ames  Research  Center 
Moffett  Field,  CA  94035-1099 


1  Air  Force  Armament  Lab 
ATTN:  AFATL/DLODL 
Eglin  AFB,  FL  32542-5000 

1  Cdr,  CRDEC,  AMCCOM 
ATTN :  SMCCR-RS P- A 

SMCCR-MU 
SMCCR  -MSI 


1  Lawrence  Livermore 

National  Laboratory 
ATTN:  Dr.  R.  Christensen 

P.O.  Box  808 
Livermore,  CA  94550 

2  Sandia  National  Laboratory 

ATTN:  Dr.  C.W.  Robinson 

Dr.  Benidetti 

Applied  Mechanics  Department 
Engineering  Mechanics  Dept. 
Engineering  Design  Division 
Livermore,  CA  94550 

1  Director 

Lawrence  Livermore 
National  Laboratory 
ATTN:  Dr.  Feng,  Mr.  Magness 
Livermore,  CA  94550 

1  Director 

Los  Alamos  National  Lab 

ATTN :  Dr .  Rabern 

Los  Alamos,  NM  87544 

2  National  Aeronautics  and 

Space  Administration 
Langley  Research  Center 
Langley  Station 
Hampton,  VA  23365 

1  Director 

National  Aeronautics  and 
Space  Administration 
Manned  Spacecraft  Center 
ATTN :  Library 
Houston,  TX  77058 

1  Director 

NASA  -  Ames  Research  Center 

ATTN ;  Tech  Lib 

Moffett  Field,  CA  94035 

1  Headquarters 

National  Aeronautics  and 
Space  Administration 
Washington,  DC  20546 
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DISTRIBUTION  LIST 


Aberdeen  Proving  Ground 

6  Dir,  USAMSAA 
ATTN:  AMXSY-D 

AMXSY-MP,  H.  Cohen 
AMXSY-G,  E.  Christman 
AMXSY-OSD,  H.  Burke 
AMXSY-TN,  D.  MetZ 

3  Cdr,  USATECOM 
ATTN:  AMSTE-CE 

AMSTE-TO-F 

13  Dir,  CRDEC,  EA 
ATTN:  SMCAR-CL 

SMCAR-CLB-PA,  M.  Miller 

SMCAR-CLB-PA,  D.  Olson 

SMCAR-CLB-PA,  J.  Huerta 

SMCAR-CLD 

SMCAR-CLY 

SMCAR-CLN 

SMCAR-CLN-D,  L.  Shaff 
SMCAR-CLN-D,  F.  Dagostin 
SMCAR-CLN-D,  C.  Hughes 
SMCAR-CLN-D,  R.  Simak 
SMCAR-CLN-D,  M.  Hayner 

2  Dir,  USACSTA 
ATTN :  Mr .  Tag 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 

This  laboratory  undertakes  a  continuing  effort  to  Improve  the  quality  of  the  reports  it  publishes.  Your 
commer.ts /answers  below  will  aid  us  in  our  efforts. 

1.  Does  this  report  satisfy  a  need?  (Comment  on  purpose^  related  project,  or  other  area  of  interest  for  which 
the  report  will  be  used.)  _ 


2.  How,  specifically,  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of  ideas, 
etc.)  _ 


3.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved, 
operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate.  _ 


4.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate  changes  to 
organization,  technical  content,  format,  etc.)  _ 


BRL  Report  Number 


Division  Symbol 


Check  here  if  desire  to  be  removed  from  distribution  list. 
Check  here  for  address  change.  _ 


Current  address:  Organization 

Address 


-FOLD  AND  TAPE  CLOSED- 


Director 

U.S.  Army  Ballistic  Research  Laboratory 
ATTN:  SLCBR-DD-T (NEI) 

Aberdeen  Proving  Ground,  MD  21005-5066 
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Director 

U.S.  Army  Ballistic  Research  Laboratory 
ATTN:  SLCBR-DD-T (NEI) 

Aberdeen  Proving  Ground,  MD  21005-9989 


